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Abstract 

The top quark may hold the key to new physics associated with the electroweak 
symmetry-breaking sector, given its large mass and enhanced coupling to the Higgs 
sector. We systematically categorize generic interactions of a new particle that couples 
to the top quark and a neutral particle, which is assumed to be heavy and stable, 
thus serving as a candidate for cold dark matter. The experimental signatures for new 
physics involving top quarks and its partners at the Large Hadron Collider (LHC) may 
be distinctive, yet challenging to disentangle. We optimize the search strategy at the 
LHC for the decay of the new particle to a top quark plus missing energy and propose 
the study of its properties, such as its spin and couplings. We find that, at 14 TeV 
with an integrated luminosity of 100 fb^ 1 , a spin-zero top partner can be observed at 
the 5<7 level for a mass of 675 GeV. A spin-zero particle can be differentiated from 
spin-1/2 and spin-1 particles at the 5<r level with a luminosity of 10 fb _1 . 



1 Introduction 



The top quark may be a window to physics beyond the Standard Model (SM). Its mass 
near the electroweak scale and its large coupling to the Higgs boson may be crucial to 
understanding the electroweak sector beyond the SM. Now that the SM-like Higgs boson 
has been observed at the Large Hadron Collider (LHC) [I] with a relatively light mass of 
about 125 GeV, the assumed "naturalness" of the Higgs sector [2] suggests the existence of a 
partner of the top quark below or near the TeV scale, motivating theories such as weak-scale 
supersymmetry, Little Higgs, and extra dimensions (either warped or universal). Vacuum 
stability of the electroweak potential also indicates the need for new physics to balance the 
large top-quark contribution. The top quark hence provides a possible early indicator of new 
physics and a good probe of a wide variety of new-physics scenarios. 

The LHC is a top factory, producing a hundred times more ti pairs from QCD processes 
than were produced at the Tevatron. Top-quark production is well understood in the SM. 
Thus any new physics contributions will be on top of a well-known and well-measured, albeit 
large, background. With the discovery era ushered in by the LHC, it would be prudent to 
keep the initial search as general as possible. 

In this work, we take a model-independent approach to searching for new physics pro- 
cesses of the form 

pp ->■ YY ->■ tiXX, 

where Y is a massive new particle with the same gauge quantum numbers as the top quark 
and X is an electrically and color neutral stable particle. The weakly interacting X could 
be a constituent of dark matter, which would manifest itself as missing energy in a col- 
lider detector. We systematically consider different spin configurations (0, 1/2, and 1) for 
the new particles Y and X. Each combination is exemplified by particles in well- motivated 
new-physics models (see the next section for details). For example, in the Minimal Supersym- 
metric Standard Model (MSSM) Y could be a scalar top and X the lightest neutralino. This 
case has been studied extensively in the literature (see, for example, Refs. [3HH]). However, 
we do not limit ourselves to specific particles in a particular model; rather, we undertake a 
general categorization, assuming merely a mass accessible at the LHC and a discrete sym- 
metry that ensure the stability of X. For simplicity, we restrict consideration to processes 
that involve only the top partner, Y, and the dark-matter candidate, X, as new particles. 

In order to distinguish experimentally between the different possibilities, one needs to 
determine the spins and couplings of the new particles Y and X. In this paper, several 
observables for this purpose are proposed and their usefulness is demonstrated in a realistic 
Monte Carlo simulation. To avoid ambiguities due to model-dependent branching fractions, 
we do not consider the total cross section in this set of variables. 

The paper is organized as follows. In section [5] we introduce the model-independent 
classification of new-physics top partners and their interactions. The production of these 
particles at the LHC is discussed in section [3j while the current bounds from collider searches 
are summarized in section HI In section [5j the expected reach of the LHC for this class of 
processes is analyzed through a detailed Monte Carlo simulation. The determination of rele- 
vant properties of the new particles, such as mass, spin and couplings, and the discrimination 
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Figure 1: Diagrams corresponding to the pair production (left) and the decay (right) of the 
color triplet Y. Double lines denote new particles, while single lines denote SM particles. 

between models are discussed in section [6j Finally, conclusions are presented in section [7J 

2 New Particles and their Couplings to the Top 

Colored particles can be copiously produced at the LHC by strong QCD interactions. Let Y 
denote a new color-triplet particle with charge +2/3. Y and its antiparticle can be produced 
at leading order in QCD by the processes shown in Fig. [T] (left). We shall not consider the 
production of a single new particle via Yukawa-type interactions: since they are strongly 
model-dependent and are subject to strong constraints from flavor physics, it is assumed 
that such vertices are forbidden by a discrete symmetry. Y decays to a new particle that is 
a color singlet, denoted X [see Fig. [1] (right)], which will show as missing energy in a collider 
experiment. 

There are four possible combinations of spins that allow a coupling between X, Y and 
the SM top quark, t. These are listed, with the relevant couplings and sample model decays, 
in Table [TJ For fermions we allow a general chirality structure. We shall henceforth refer to 
these scenarios as models i, ii, iii, and iv. 

Let us elaborate on the unusual case in which Y is a vector color triplet, possibly arising 
as a bound state from strong dynamics or from a special kind of supersymmetric model [12j. 
The kinetic term is 

C*in = ~(F^F^, F IM/ = D^Y v -D v Y li , (1) 
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Table 1: Quantum numbers and couplings of the new particles X and Y, which interact 
with the SM top quark, t. In the last column, t and Xi are the scalar top and lightest 
neutralino in the MSSM, respectively fWj . t KK , ^kk? and 7_h" i kk are the first-level Kaluza- 
Klein excitations of the top, the photon, and an extra-dimensional component of a photon, 
respectively, in universal extra dimensions (UED) [ilf . Finally, Q is the vector superpartner 
in a supersymmetric model with an extended gauge sector fl^ . 



where _D» = du — igT a G a . Then the y-y-gluon interaction term is 



£yyg = ^g(T a h ((cPy/ - d%n(G*Y ui - G a u Y^ - h.c.) . 



(2) 



The resulting Feynman rule is 
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Likewise, the y-y-gluon-gluon interaction term is 

9 2 
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(4) 



The resulting Feynman rule is 




-ig 2 {{T c T d + T d T c ) k ,g TX g^ - (T c T d ) k ^g x " - (T d T c ) k ^g x °) 



(5) 
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Figure 2: Production cross sections for pp — > YY at the LHC for 8 TeV (left) and 14 TeV 
(right), as a function of the mass my, for different Y spins. The leading SM backgrounds 
are indicated by horizontal lines. 



3 Color- Triplet Top-Partner Production 

The dominant modes for production of the top partner in hadronic collisions are the QCD 
subprocesses 

qq,gg^YY. (6) 

We restrict ourselves to the first- and second-generation quarks q = u,d,c,s and use the 
CTEQ 6L1 parton distribution functions (PDFs) 1 1 3 j . with the factorization scale set to 
my. For my ~ 200-1000 GeV, the dominant subprocess is gg — > YY, which is about one 
order of magnitude larger than uu, dd — > YY. The channels cc, ss — > YY are suppressed by 
roughly one additional order of magnitude. 

The total QCD production cross section at the LHC as a function of the mass of the Y 
is shown in Fig. [21 for the cases in which Y has spin 0, 1/2, and 1. The plots include next- 
to-leading order (NLO) and resummed next-to-leading logarithmic (NLL) QCD corrections 
for the scalar Y [H] and NLO and NNLL corrections for the fermionic Y [15J. The QCD 
corrections for vector Y production have not yet been calculated; we use the i^-factor for 
the scalar Y (1.77 at ^/s = 8 TeV, 1.56 at ^/s = 14 TeV), since the two cases share similar 
leading-order Feynman diagrams. 

The cross section for the fermion is about an order of magnitude larger than that of the 
scalar, because of the fermion's extra spin degrees of freedom and threshold effects. In the 
s-channel, the scalar is produced as a p-wave with a velocity dependence of o ~ j3 3 , whereas 
the fermion is produced as an s-wave with a ~ /3. Thus, the ratio of the cross sections of 
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Table 2: Experimental bounds on the mass of particle Y for different spins, Jy, under the 
assumption my ^> mx ■ These estimates are based on the ATLAS results from Refs. JI9J . 

the fermion and scalar is larger at small values of /3, that is, when the mass of the Y is large. 
This relative enhancement of the fermionic Y Y production is particularly pronounced when 
the Y particles are produced mostly near the threshold limit. Note that, although the curves 
for the vector and scalar appear to be parallel on the logarithmic scale, their ratio varies 
from about 34 to 16 in the mass range shown. 

4 Current Bounds from the Tevatron and LHC 

As Fig. [1] indicates, the top-quark partner, Y, decays to a top quark plus a neutral particle, 
X. The discrete symmetry implies that X is stable and leads to missing-energy events. 
Thus, the signal is tt plus missing energy. Searches for supersymmetric scalar tops at the 
Tevatron [H1[T7] and the LHC [TSH2T)] put constraints on the allowed parameter space for 
the class of processes considered here. Additional, though generally weaker, bounds also 
arise from general searches for signals with jets and missing energy \21\ . 

Currently, the strongest constraints arise from scalar top searches at ATLAS using 
4.7 fb _1 of data taken at y'i = 7 TeV [19J. For my 3> mx, they put a lower bound 
my > 500 GeV on a scalar Y. By taking into account the different production cross sections 
for scalars, fermions, and vectors (see Fig. [2]), one can translate the results of Refs. [19J into 
a limit of my > 650 GeV for a fermionic top partner and my 730 GeV for a vector top. 

The bounds are summarized in Table [2j It should be pointed out that the limits for 
fermionic and vector Y are simple estimates from theoretical considerations. For more 
robust results, a detailed experimental analysis of these scenarios needs to be performed. 

For larger values of mx, that is, smaller mass differences my — mx, the limits become 
weaker. The excluded region in the my-mx mass plane for scalar Y particles will be shown 
in the next section (see Fig. HJ) . 

5 Signal Observability at the LHC 

As the previous section discusses, we consider new physics signals of the type tt + If}. For 
the leading channel, in which the top quarks decay hadronically [3], the signal receives large 
backgrounds from SM processes with multiple QCD jets. To suppress QCD backgrounds, 
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we consider the semileptonic channel [I], in which one of the tops decays hadronically and 
the other decays leptonically, namely, 



pp -> YY -> £X tX -> 6jij 2 6TP* XX + h.c. (£ = e, /i). (7) 

This channel is beneficial because of its sizeable branching fraction and the identification of 
both t and i. The dominant background processes are 

tt, ttZ (with Z ->■ i/P), and j (with W ->• . (8) 

The cross sections for the first two backgrounds (without branching fractions) are shown in 
Fig. [2] as horizontal lines, including NLO corrections for tiZ [22] and NLO+NNLL effects 
for tt |23]. 

The separation of signal and background in the semileptonic channel has been studied 
previously in the literature [Hl5|[T6 t[T8l420] . Here, we reanalyze the signal selection with the 
purpose of developing optimized selection cuts in a phenomenologically realistic simulation 
setup. Our signal selection follows the strategy of Ref. [I], but we include QCD parton 
showering and detector smearing effects. As a result, we find that we need to adjust the 
choice of cuts to account for the effect of QCD radiatiorfl 

Jets have been clustered via a cone algorithm with cone size 0.4. To simulate detector 
resolution effects, we have smeared the jet energy with a Gaussian distribution of width 
0.5 x \/E, where E is the jet energy in units of GeV. A 6-tagging efficiency of 70% [24J has 
been assumed. We have applied the following set of cuts to identify the signal signature and 
reduce the SM backgrounds. 

exactly one lepton £ — e, fj, with > 20 GeV, < 2.5; 
at least two light jets with E J T > 25 GeV, \rjj\ < 2.5; 
. , exactly two b-tagged jets with E\ > 30 GeV, |7#,| < 2.5; , . 

U ARjj, AR bj , AR bb > 0.4, AR ej = AR ib = 0.3; 

70 GeV < m j:i < 90 GeV, 120 GeV < m^ ad < 180 GeV; 
E > 25 GeV. 

Here b and j stand for a jet with or without a 6-tag, and E l T and rji are the transverse energy 
and pseudorapidity of object i. AR = a/ (Ar/) 2 + (A0) 2 describes the angular separation 
between two jets. m^ ad is computed from either the bjj or the bjj invariant mass, namely, 
whichever yields the value closer to the true top-quark mass, m t , in a given event. ]fi is the 
missing transverse energy. 

Events for the partonic signal process and ttZ background have been generated with 
CalcHEP 3.2.5 [25] and passed to Pythia 6.4 [26] for parton showering and jet clustering. 
The tt background has been simulated with PYTHIA. It was shown in Ref. [1] that the Wbbjj 



* Very recently, several papers have appeared that pursue a similar goal in the context of the MSSM, using 
traditional selection cuts [BJ[5] and top-jet tagging techniques [7] . Our results for the signal observability are 
comparable to Refs. [HE], but significantly better than Ref. [8]. 
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Figure 3: Differential cross sections for the SM backgrounds ti (red, solid) and tiZ (blue, 
dashed) and the signal scalar Y with (My,Mx) = (600, 10) GeV (cyan, dash-dotted) before 
and after the cut M^' miss > 90 GeV. Distributions after the M^' miss cut are shown in bold 
lines. 

background can be reduced effectively with invariant-mass cuts on the jj for a W selection 
and bjj for a top-quark selection. We have thus neglected this process in our simulation. 

With the set of cuts in ([9]), which we shall refer to as Cut( 1), a. good signal-to-background 
ratio is achieved for small values of my, when the YY production cross section is large. For 
larger values of my, additional cuts are required to suppress the SM background sufficiently. 
It turns out that the following two variables are useful for this purpose: the missing transverse 
energy, ]fi, and the transverse mass of the lepton-missing- momentum system, 

i\4 miSS = V /(E ff + ^-(p, T + F i T )2. (10) 

The optimal cut values depend on the collider energy: 

• For y/s = 14 TeV, the choice 

Cut(2h) = Cut(l) plus ^> 350 GeV and M^' miss > 90 GeV 

has been found to be effective for my ~ 600 GeV. It can be understood as follows. 
The cut il4' miss > 90 GeV is necessary because a large amount of missing energy in 
the SM backgrounds corresponds to neutrinos from the leptonic decay of the W boson. 
From Fig. [3] one can see that the M^' miss cut reduces the SM backgrounds dramatically, 
especially in the region. However, the signal events remain virtually the same 

after this cut: only those in the low- -$5 region are slightly affected. Moreover, for the 
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Figure 4: Expected statistical significance for model i, as a function of the masses of X 
and Y . The left panel corresponds to ^/s = 8 TeV and C = 20 fhr 1 , while the right panel 
corresponds to \fs = 14 TeV and C = 100 fb" 1 . The dashed line shows the current exclusion 
limit at the 95% confidence level from Ref. \19]j . 

signal there is a plateau between 200 and 400 GeV in Fig [3j The cut ft > 350 GeV 
has been chosen because above 350 GeV the backgrounds are suppressed considerably. 
In practice, we have applied either Cut(l) or Cut(2h), whichever produces the larger 
statistical significance S/yfB for a given parameter point (my, mx)- Here, S and B 
denote the number of signal and background events after cuts. 

• For yfs = 8 TeV, we have used either Cut( 1 ), or 

Cut(2l) = Cut(l) plus E> 200 GeV and M^' miss > 145 GeV, 

or 

Cut(2V) = Cut(l) plus ^> 300 GeV and M^' miss > 185 GeV, 

whichever results in the largest significance. Cut (21) and Cut (21') have been optimized 
for my ~ 400 GeV and my ~ 500 GeV, respectively. 

Because the signal cross section is lower for yfs = 8 TeV than for y/s = 14 TeV, we 
have lowered the ft cut to ensure that a sufficient number of signal events passes. 
However, the looser ft cut also results in a larger background event yield, so that it is 
advantageous to apply a stronger cut on M^ miss to improve the signal significance. 

Figure H] shows the statistical significance that can be achieved with these cuts for model i 
(scalar Y and fermionic X), for different values of my and mx- The significance is deter- 
mined according to S/ \/B if B > 10, whereas Poisson statistics is used for very low event 
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spin-0 


spin- 1/2 


8 TeV 
14 TeV 


480 GeV 
675 GeV 


660 GeV 
945 GeV 



Table 3: The 5a discovery reach for spin-0 and spin- 1/2 top partners at 8 and 14 TeV with 
integrated luminosities of 20 and 100 fb~ l , respectively, mx = 100 GeV is assumed. 

yields (B < 10). For the other model combinations, ii-iv, the statistical significance can be 
obtained by scaling the values in Fig. H] with the production cross sections in Fig. |2j 

As the figure shows, the statistical significance is relatively large in the following two 
regions of the mass plane: 

1. Small values of my, in which case the signal selection efficiency is almost independent of 
mx- Here, the YY pair, which recoils against an initial-state jet, is typically produced 
with a sizeable boost. This boost leads to a fairly large missing momentum, which 
helps to discriminate the signal from the ti background. 

2. Moderately large values of my, my ^ 600 GeV, and small values of mx- For these 
values of my, the YY pair is produced mostly at rest and the signal selection becomes 
difficult for small mass differences my — mx, when the top quark from the decay 
Y — > tX is quite soft. 

In comparison with Ref. [4], we obtain somewhat lower values for the significance S/y/B, 
as a consequence of having performed a more realistic simulation that includes QCD radiation 
(through parton showering) and jet smearing. These effects make it more difficult to devise 
clean kinematic selection variables for the signal and result in more background from the tail 
of smeared distributions. We have also explored the mass reconstruction scheme proposed 
in Ref. [4] and the variable Mti [2E]- We have found them to be useful in certain respects 
and complementary to the combination of our cuts. Further optimization would depend on 
detailed (experimental) simulations, which we leave for future studies. 

In summary, we have found that, at 14 TeV with an integrated luminosity of 100 fb _1 , 
a scalar top partner can be observed at the 5a level (or better) for a mass up to 675 GeV 
if Mx = 100 GeV. This translates into 945 GeV for a spin-1/2 top partner. At 8 TeV with 
an integrated luminosity of 20 fb _1 , it is possible to achieve a 5a discovery for a scalar top 
with a mass up to 480 GeV. This corresponds to 660 GeV for a spin-1/2 top partner. These 
results are summarized in Table [3j 

6 Determination of Model Properties 

6.1 Masses 

The independent determination of the Y and X masses in processes of the type 

pp^YY -> ffXX, (11) 
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where / is a SM fermion, is a difficult problem because of the lack of kinematic features for 
the under-constrained system. Several methods have been proposed in the literature [2ZH2I], 
either based on global event variables such as M eff = J2 i€vis Vt,% + $t, on the variable 

M T2 = min PT Xi+PT X2= ^ T jmax (M,f > X \ M^~' X2 ) } [28] or variants thereof, or on likelihood 

fits to the complete event information [2S]. It was found that, for my ~ (9(300 GeV) and 
a sample of a few tens of thousands of signal events at y/s = 14 TeV, the mass difference 
my — rrix can be determined to a precision of a few per cent, while the absolute mass scale 
has an uncertainty of roughly 20-30% [29l[30]. If the YY state could arise from the decay 
of a new resonance of known mass, it would help to constrain the kinematics and thus to 
determine the masses of Y and X as well [3T] • More details can be found in the cited papers. 



6.2 Spin 

The spin of the Y particle can be probed through the characteristics of the YY production 
process. For instance, the YY production cross section strongly depends on the spin [32J. 
However, unknown model-dependent branching fractions and the mass uncertainty of or- 
der 30% can lead to ambiguities in the determination of the spin from the measured total 
production rate. Instead, one can largely avoid such problems by investigating the shape 
of suitable differential distributions. In particular, the two variables described below are 
effective for this purpose. 

(1) Scalar and fermion Y pair production can be distinguished with the observable 

tanh(Ay tf /2), Ay e = \y bjj - y bi \, (12) 

which is constructed from the rapidities of the visible decay products of the hadronically 
decaying and the leptonically decaying top quarks. In general, there is a combinatorial 
ambiguity in identifying the 6-jets and light-flavor jets as the decay products from one 
of the two top quarks. Given our event reconstruction scheme discussed in the previous 
section, one can resolve this ambiguity by assuming that the hadronically decaying top 
quark is made up from the two light-quark jets and the 6-jet for which m^jj is closest 
to m t . The remaining 6-jet and the lepton are then identified as the decay products of 
the other top quark. 

The variable in (|12p is closely related to the proposal by Barr in Ref. [33], tanh(|m- — 
r}f\/2), where r)f is the pseudorapidity of the SM fermion from the decay Y —> Xf. 
This variable approximately traces the production angle 9* between one Y and the 
beam axis in the center-of-mass frame. In the qq — > YY channel, the 9* distribution 
has a clear dependence on the Y spin, as can be seen in the formulae 

-[qq — > YY] oc 1 — cos 2 ^*, for scalar Y (spin 0), (13) 



d cos 9* 

, —^ [qq -> YY] oc 2 + /3 y (cos 2 9* - 1), for fermionic Y (spin |), (14) 

& COS v 
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Tan h (Ay/2) 

Figure 5: Distribution of tanh(A?/ t f/2) for the different top-partner scenarios listed in 
Table for my = 300 GeV, mx = 100 GeV, and \fs = 14 TeV. For comparison, all 
distributions have been normalized to unity. 



where (3y is the velocity of the produced Y particles. The difference stems from the fact 
that scalars are produced in a p-wave, whereas for fermions the s-wave contribution is 
dominant. In contrast to Ref. [33], the definition (fl~2|) is based on the rapidities rather 
than the pseudorapidities, to account for the fact that the produced top quarks are 
massive. 

In the physical process pp — > YY, only a subdominant fraction of the events originates 
from qq annihilation, but as can be seen in Fig. |5]the effect is still noticeable (compare 
the lines for model i with the other cases). 

(2) There is no appreciable difference between fermionic and vector Y pair production 
in the tanh(Ay t f/2) distribution. However, these two cases can be disentangled by 
means of a variable that measures the effective hard scattering energy [31]. One such 
observable is the effective mass, a scalar sum over momenta: 

M eS = Pta + (15) 

igvis. 

where the sum runs over all visible objects (jets and leptons in this case). 

The usefulness of this variable follows from the fact that the partonic cross section for 
the pair production of massive vector particles grows with the partonic center-of-mass 
energy like s / m Y , whereas for fermions it has the usual 1/s behavior in the high-energy 
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Figure 6: Distribution of M e Q for the different top-partner scenarios listed in Table HJ for 
my = 300 GeV, mx = 100 GeV, and y/s = 14 TeV. For comparison, all distributions have 
been normalized to unity. 



limit. In fact, for very large values of s the vector YY production would violate the 
perturbative unitarity limit, a sign that additional new massive resonances will appear 
and modify the production amplitude. The presence of such resonances is generally 
expected in dynamical models such as the one proposed in Ref. [12J. However, their 
masses may be beyond the reach of the LHC, depending on the value of s when the 
unitarity limit is eventually reached. We have estimated this limit for vector Y pair 
production, conservatively assuming s-channel dominance in obtaining s > (4.8 TeV) 2 . 
Thus, one can assume that the new heavy resonances have masses of (9(5 TeV), which 
make their potential contribution to the process pp — > YY completely negligible, since 
only a fraction of order 10~ 6 of events have partonic center-of-mass energy of this size 
or larger at =14 TeV. 

The M e fj distribution is shown in Fig. [6] for the different model scenarios. A very 
distinctive difference can be observed between the cases of vector and fermionic or 
scalar top partners. 

Let us quantify the discriminative power of these variable in an example: the process 
pp YY ->■ tiXX for m Y = 300 GeV and m x = 100 GeV at ^ = 14 TeV. The simulation 
and event selection have been performed as described in the previous section, with Cut(l) 
in Eq. (J9]). We have assumed the cross section for scalar Y pair production (model i). Note 
that this scenario is not ruled out by current LHC results. We have not used the total event 
rate for model discrimination, to avoid ambiguities due to unknown branching fractions. 
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We have carried out the discrimination between two different spin assignments by com- 
puting the x 2 value for the binned tanh(Ay t f/2) and M eS distributions, using three bins in 
both cases. The result can be expressed in terms of the integrated luminosity £ 5cr necessary 
for achieving a 5cr statistical significance: 



14 TeV: scalar Y versus fermion Y: 




= 9.4 fb~ 


-i 


scalar Y versus vector Y: 




= 0.8 fb" 


-i 


fermion Y versus vector Y: 


£>5a 


= 0.7 fb" 


-i 



For the current 8 TeV run, a 5a discrimination requires the following integrated luminosities: 

8 TeV: scalar Y versus fermion Y: £ 5(T = 72 fb -1 , 
scalar Y versus vector Y: C 5a = 8.1 fb -1 , 
fermion Y versus vector Y: £ 5(T = 5.2 fb -1 . 

The numbers refer to the purely statistical significance. However, at this level of precision, 
systematic errors may be important. A potentially large systematic effect stems from the 
uncertainty of the new-particle masses, my and mx- While the mass difference my — mx 
can be determined rather precisely, the overall mass scale can be measured with only 20-30% 
accuracy; see section 16.11 We have estimated the effect of this uncertainty by comparing 
two event samples with (my, mx) = (300, 100) GeV and (my, mx) = (400, 200) GeV, which 
differ in my by roughly 30%. We have found that this mass uncertainty reduces the statistical 
significance of the spin discrimination by about 20%; the values of £ 5(T that account for this 
systematic error are about 50% greater than those quoted above. 

In conclusion, the determination of the spin of the top partner, Y, is possible with very 
moderate amounts of data. On the other hand, the distinction between models ii and iii, 
which both have a fermionic Y but differ in the spin of the singlet X, is much more difficult. 
After surveying more than a dozen different kinematic variables based on the top-quark 
momenta, we found no significant difference between scenarios ii and iii for any of them. 
This finding agrees with the results of Ref . [34] . 

However, more information can be obtained from observables that are sensitive to the 
top-quark polarization, as will be discussed next. 

6.3 XY Couplings 

The chirality structure of the decay Y — > tX (that is, the relative contributions of left- and 
right-handed chiral couplings) leaves an imprint on the polarization of the top quark, which 
can be analyzed through angular distributions of the top-quark decay products. This method 
is particularly effective when the mass difference between Y and X is large (my — m x 3> mt), 
so that the top quark is energetic and therefore the helicity is preserved, reflecting the 
chirality. For instance, one can look at the angle 8' b (8' £ ) of the b quark (lepton) with respect 
to the top-quark boost direction in the top rest frame. Because the b quark is always left- 
handed, it is emitted predominantly in the forward direction (cos 9' b > 0) if the top quark 
is left-handed, but mostly in the backward direction (cos 6^ < 0) if the top quark is right- 
handed. 
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model ii 




-0.5 0.0 0.5 

cos(e' ( ) = -cos(^) 



0.5 0.0 0.5 

cos(6>;) = -cos(^) 



1.0 



Figure 7: Parton-level angular distribution of the top-quark decay products in the b£ rest 
frame, for the decay chain Y — > Xt — > Xbi^v. The four panels show the results for the four 
scenarios in Tabled for the two coupling choice — 1, clr — (black) and a L = 0, a R = 1 
(red). The input mass parameters are my = 400 GeV and mx = 10 GeV. The distributions 
have been normalized to unity. 



In practice, even if the top quark is produced fully polarized in the decay Y — > tX, some 
of the polarization is washed out by the mass of the top, but the cos 6' b < distribution will 
still exhibit a characteristic difference between left- and right-handed XYt couplings. 

In the following, we shall illustrate this behavior using a parton-level simulation with 
CalcHEpQ. We shall focus on the leptonically decaying top quark, since it has a cleaner 
final state. The top-quark rest frame cannot be reconstructed because of the unobserved 
neutrino momentum, so we analyze the angular distribution in the rest frame of the visible 
bl system instead. The results are shown in Fig. [7J 

As Fig. [7J shows, the distribution is skewed to smaller values of cos 6' t or, equivalently, 
larger values of cos9' b in the case of a left-handed XYt coupling (black curves) than in 
the right-handed case (red curves). For a mixed case with non-zero left- and right-handed 
components, one obtains a distribution that lies between the black and red curves. This 
qualitative behavior is the same for all four spin combinations in Table [TJ although they 



tA more realistic simulation at the level of the previous section, including parton showering and signal 
selection cuts, would require the modification of Pythia to include top-quark spin correlation effects, which 
we have not attempted to carry out. 
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differ from each other in the detailed shape of the distribution. In particular, cases ii and 
iii have distinctly different shapes; hence, the analysis of this observable may allow one to 
determine not only the chirality of the XYt coupling but also the spin of the X particle. 
Such a determination is not possible with observables that treat the top quarks as basic 
objects. 

Furthermore, one can probe the chirality even with limited statistics by using two bins 
and forming the asymmetry 

Af„\ - ff ( cos ^ > x ) ~ ^(cosfl; < x) 
a[cosUf > x) + a{cosu e < x) 

From Fig. [7J one can see that when x is about —0.5 A(x) will be most sensitive to the 
chirality of the coupling. Table H] shows the asymmetry A(— 0.5) for models i-iv with two 
choices of the masses my and mx- The usefulness of A{— 0.5) in determining the coupling 
is enhanced by its relative insensitivity to the spin and mass combinations. 





ITLy = 


400 GeV, m x = 


10 GeV 


my = 


300 GeV, 


m x = 


100 GeV 








Model 






Model 




a>L, 


a>R 


i 


ii 


iii 


iv 


i 


ii 


iii 


iv 


1, 





-0.10 


0.02 


-0.10 


-0.03 





0.15 


0.04 


0.10 


o, 


1 


0.68 


0.55 


0.68 


0.61 


0.54 


0.39 


0.50 


0.45 


1, 


1 


0.29 


0.28 


0.29 


0.29 


0.28 


0.27 


0.28 


0.27 



Table 4: Asymmetry A{— 0.5) for models i-iv and two choices of the masses my and mx- 



7 Conclusions 

The SM-like Higgs boson has been observed at the LHC with a relatively light mass of about 
125 GeV. The "naturalness" argument of the Higgs sector suggests the existence of partners 
of the SM particles, especially the heavy top quark. The top quark may thus hold the key 
to new physics associated with the electroweak symmetry-breaking sector, because of its 
enhanced coupling to the Higgs sector. In this paper, we have systematically categorized 
the generic interactions of a new particle that couples to the top quark and a stable neutral 
particle, which serves as a candidate for cold dark matter. We have considered all possible 
assignments of spin 0, | and 1 for either of the two new particles. 

In the search for new physics involving top quarks and its partners at the LHC, the 
experimental signatures may be distinctive, but challenging to disentangle. Pair production 
of the massive top partners leads to a signature of a tt pair plus missing energy, which is 
difficult to separate from the large SM tt background. We have presented a set of optimized 
selection cuts for isolating this new physics signal at the 8 and 14 TeV runs of the LHC. 
We have found that, at 14 TeV with an integrated luminosity of 100 fb _1 , a spin-zero top 
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partner can be observed at the 5a level up to a mass of 675 GeV, while for a spin-| top 
partner the reach extends to 945 GeV. 

If a process of this type is discovered at the LHC, it will be imperative to determine 
the spins and couplings of the new particles, in order to understand the underlying physics 
mechanism. We have proposed a strategy to extract these properties from experimental 
data by means of suitable differential distributions of the final-state products. With this 
approach, a spin-0 top partner with mass of about 300 GeV can be discriminated from spin- 
| and spin-1 particles at the 5a level with a luminosity of 10 fb _1 at 14 TeV. Furthermore, 
the structure of the coupling that mediates the decay of the top partner into a top quark 
and a massive neutral particle can be analyzed by measurement of the polarization of the 
final-state top quarks. This method allows one to distinguish clearly between left-handed, 
right-handed, and vector couplings. Most importantly, the proposed observables for spin 
and coupling determination are insensitive to unknown branching fractions and depend only 
mildly on the masses of the new particles. 

In conclusion, the LHC will allow us to observe and study top partners with mass up to 
about 1 TeV. This program will shed light on the interplay of the Higgs-boson and top-quark 
sectors and may elucidate the concept of naturalness. 
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